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ABSTRACT 


A  mathematical  model  is  developed  for  determining  an  optimum 
submarine  tender  load  list  for  the  POLARIS  (FBM)  weapon  system. 
Military  essentiality  of  spare  parts  and  the  stowage  capacity  of  the 
tender  play  prominent  roles  in  the  development  of  the  model.  The 
measure  of  effectiveness  which  is  maximized  is  the  average  number 
of  demands  filled  over  a  given  period  of  time.  This  is  equivalent  to 
minimizing  the  number  of  stockouts  of  an  item.  An  example  of  the 
use  of  the  model  is  included  for  the  case  where  the  demand  for 
items  is  characterized  by  the  normal  distribution.  The  paper  in¬ 
cludes  a  computer  flow  chart  and  program  to  aid  in  implementing 
the  model. 

The  authors  wish  to  express  their  appreciation  for  the 
assistance  and  encouragement  given  them  by  Associate  Professor 
Julius  H.  Gandelman  of  the  U.  S.  Naval  Postgraduate  School. 

The  authors  are  also  indebted  to  Miss  Patricia  Hoang  for  her  assis¬ 
tance  in  programming  the  model  on  the  CDC  1604  computer. 
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GLOSSARY  OF  TERMS 


zi 

xi 

*i(xi) 

eM 

b(zj) 

s(zj) 

X 

ci 

c 

c 

Mi 

mi 

P 


a' given  item  where  i  =1,2,3 . N  (the  entire 

list  of  spare  parts  is  drawn  up  and  then  a  number 
is  assigned  to  each  spare  part  to  identify  it). 

amount  of  item  i  to  be  stocked. 

tVi 

demand  for  the  il  item. 

probability  that  an  amount  Xj  will  be  demanded 
during  the  period  under  consideration. 

expected  demand  for  item  i  ;  =  mi. 

expected  shortage  of  item  i. 

th 

average  number  of  1  n  item  supplied  to  meet  demand. 

Lagragian  multiplier. 

cube  (volume)  of  the  i ^  item. 

total  cube  available  for  spare  parts  stowage  on  tender, 
total  cube  associated  with  any  A  for  which 

th 

essentiality  (military  worth)  number  assigned  to  i 
item. 

th 

deviation  of  i  item  from  the  mean, 
mean  of  the  i^  item. 

th 

characteristic  number  associated  with  the  i  item. 


v 


1.  Introduction. 


The  trend  toward  the  development  and  utilization  of  highly  com¬ 
plex,  sophisticated  weapon  systems  within  the  Navy  has  generated 
an  ever-increasing  requirement  for  more  efficient  logistical  support. 
This  paper  is  an  attempt  to  develop  a  mathematical  model  that  can 
be  used  as  a  decision  aid  for  one  aspect  of  the  logistics  problem 
associated  with  the  POLARIS  system.  The  model  provides  a  method 
for  determining  the  items  and  associated  quantities  that  should  be 
stocked  aboard  the  supporting  tender. 

The  requirement  for  more  effective  tender  load  lists  is  even 
more  evident  when  viewed  in  conjunction  with  the  POLARIS  mobile 
base  concept. 

In  the  event  of  enemy  attack,  the  Navy's  immediate  POLARIS 
retaliatory  capability  will  have  been  dependent  on  the  ability  of  the 
supporting  tender  to  provide  the  necessary  SSB(N)  replenishment 
and  repair  requirements.  This  aspect  of  the  tender  function  is  fur¬ 
ther  complicated  by  the  limitations  of  stowage  space  available  for 
stocking  the  required  spare  parts.  As  a  result  of  this  space 
limitation,  tender  load  lists  must  reflect  SSB(N)  requirements,  and 
tender  re- supply  procedures  must  be  closely  coordinated  and 
evaluated  to  insure  that  SSB(N)  requirements  are  being  met  in  an 
effective  and  efficient  manner. 
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The  concept  of  military  worth  is  incorporated  in  the  model  to 
reflect  the  assignment  of  priorities  between  the  various  items. 

These  essentiality  numbers  are  an  attempt  to  define  the  military 
worth  of  a  particular  item  with  reference  to  the  performance  of  the 
integrated  weapon  system.  Essentiality  factors  will  be  assigned 
arbitrarily;  however,  it  should  be  apparent  that  actual  MEC  numbers 
could  have  been  utilized  for  illustrative  purposes. 

The  measure  of  effectiveness  to  be  maximized  will  be  the  average 
number  of  demands  for  spare  parts  that  are  filled  by  the  tender  over 
some  period,  as  a  function  of  the  military  essentiality  of  the  various 
items.  This  is  equivalent  to  minimizing  the  number  of  stockouts 
of  an  item. 

It  is  obvious  that  if  the  demand  for  an  item  was  known  exactly, 
the  problem  would  be  simplified  considerably  subject  to  the  space 
limitation  of  the  tender.  Although  demand  cannot  be  predicted  with 
certainty  it  may  be  possible  to  estimate  the  probability  distribution 
of  demand  for  each  item  on  the  basis  of  past  usage  data.  In  the  event 
that  usage  data  are  not  available,  the  model  can  be  utilized  by  as¬ 
suming  a  probability  distribution  of  demand  for  the  item(s),  based 
upon  experience  and  judgment.  This  would  be  the  situation  encoun¬ 
tered  in  developing  an  initial  tender  load  list  for  a  new  class  of 
submarines.  As  usage  data  become  available  the  tender  load  list 
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can  be  continually  updated  to  reflect  actual  demand,  utilizing 
automatic  data  processing  equipment.  Appendix  III  contains  a  com¬ 
puter  flow  diagram  and  program  that  will  accomplish  this  objective. 

Although  this  paper  is  primarily  concerned  with  presenting  a 
method  for  determining  an  optimal  tender  load  list,  certain  aspects 
of  the  tender  re- supply  problem  are  included  for  consideration. 
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2. 


Assumptions 


The  following  simplifying  assumptions  are  introduced  in  order 
to  illustrate  the  proposed  method. 

1)  The  model  is  concerned  with  spare  parts  and  does  not 
consider  dry  stores  and  other  consumable  items.  As  a  result,  C 
refers  to  the  cube  available  for  spare  parts  stowage. 

2)  is  a  continuous  variable,  and  its  associated  probability 
density  function  f^(x^)  possesses  the  properties  required  for  the 
existence  of  derivatives  and  integrals. 

3)  The  demand  for  each  item  is  independent  of  the  demand 
for  every  other  item. 

4)  The  probability  distribution  of  demand  is  known  for  each 
item  over  the  time  period  considered. 

5)  The  list  of  items  considered  reflects  only  those  for  which 

demand  has  been  generated.  If  demand  for  an  item  stops  ,  then  the 
item  is  removed  from  stock  and  a  new  load  list  is  generated.  Con¬ 
versely,  if  demand  for  a  new  item  is  generated,  it  is  added  to  the 
load  list  and  a  revised  Z  =  |z^  z2)'  •  ZN  |  determined  by  means 

of  the  model. 

6)  The  development  of  the  model  is  independent  of  the  re¬ 
supply  phase  of  the  tender  operation.  The  load  list  determined  is  a 
function  of  demand  over  some  period  during  which  the  tender  must 
provide  support  without  being  re- supplied.  Theoretically  the  tender 
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is  completely  depleted  of  spare  parts  at  the  end  of  the  period  con¬ 
sidered.  The  effect  of  re- supply  as  a  function  of  the  desired  supply 
objective  will  be  considered  in  a  later  section. 

7)  The  utility  values  defined  in  Appendix  I  are  additive. 
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3.  Presentation  of  the  Model. 

This  section  presents  the  essential  concepts  of  the  proposed 
model.  Appendix  I  contains  a  detailed  mathematical  formulation  of 
the  results  of  this  section. 

The  expectation,  ECU,  of  x.  with  respect  to  the 

probability  function  PL.*3  gives  the  following  expression  for  the 
average  amount  of  item  i  demanded: 

(!)  EC*il  =  fi  (*i)  dxi» 

where  Xi=  demand  for  the  i^1  item. 

(The  limits  of  integration  range  from  0  to  OO  since  only  items 
for  which  a  demand  is  generated  are  considered) 

Consider  the  following  situations: 


(a)  x{  <  zi  , 

(b)  xi  >  z^, 


where  zj  =  amount  of  item  i  to  be  stocked. 

If  demand  for  an  item  is  less  than  the  amount  of  the  item  that 
isstocked  then  the  tender  can  meet  demand.  However,  if  demand 
exceeds  the  stock  level  then  only  an  amount  zi  can  be  supplied.  If 
this  condition  occurs  (x^  -  z^)  represents  the  shortage  (backorder) 

associated  with  that  item.  The  expected  shortage  of  an  item  can  be 
expressed  as  a  function  of  the  amount  of  the  item  stocked. 

(2)  b^z^  =  J  (xi— z^  fj(xi)  dx^  For  x^zp 


0 


otherwise  . 


The  expected  value  of  the  amount  actually  supplied  to  meet 

demand  is  given  as  a  function  of  the  amount  stocked  on  the  tender. 
This  results  in  the  following  expression: 

(3)  Si£Zi3  =  J*  J  34  fi(x.)  dx{  +  Zi  fiUiJdxi 
Manipulation  of  equations  (1-3)  confirms  the  intuitive  concept 

that  the  expected  demand  must  equal  the  expected  demands  filled 

plus  the  expected  shortages. * 


(4)  ECxil  =  +  bi  ^3 

As  mentioned  previously,  the  measure  of  effectiveness  is  a 
function  of  the  military  essentiality  of  the  item  as  well  as  the  de¬ 
mands  filled  by  the  tender.  This  requirement  is  met  by  defining 

N 

M(z, ,  .  .  zn)=  2M.  s^(z-)  in  accordance  with  assumption  (7). 

1  c*  1 

The  problem  is  now  one  of  determining  the  set  ^=J|  •*'2.^ 

such  that 


is  maximized  subject  to  the  constraint  cf  c-  z-  =  C  ,  z^>0. 


Equation  (5)(below)  is  derived  in  Appendix  I  and  is  a  general 
result  of  the  considerations  mentioned  in  this  section.  This  equation 
provides  the  means  for  determining  an  optimum  tender  load  list. 

(5)  f  fi(xi)  dxi  =l-\(2L),  i=  1,  2 . N, 

*0 

where  A  is  a  lagrangian  multiplier  with 

M  * 

values  in  the  range  min  (  — —  )  >  A  —  O 
subject  to:  ^ 

_  Jr  =  C  n  ^i>o. 

^Appendix  I 
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4.  General  procedure  for  utilizing  the  model. 

The  general  procedure  for  solving  equation  (5)  is  to  assume  a 

X*A  ' 

corresponding  to  min  (  ) ,  and  compute  the  set  Z 

by  using  the  appropriate  cumulative  distribution  function  table.  The 
C1  obtained  (where  C'  results  from  using  a  A  other  than  that  for 
which  when  the  are  substituted  into  the  constraint 

equation  (  21  -  C-  )  determines  the  minimum  required  stowage 

space  aboard  the  tender  as  a  function  of  demand  over  the  demand  period 
associated  with  f^(x-).  (As  X  increases,  C1  decreases.  Since 
min  (  )  is  the  largest  permissible  value  of  X  >  the  corresponding 

C/is  a  minimum,)  It  should  be  noted  that  for  any  \  in  the  range 
min  X  —  O  ,  every  ^  0.  This  automatically  sets  a  mini¬ 

mum  stowage  requirement  on  the  tender  of  at  least  the  sum  of  the 
cube  of  each  item. 

Under  the  assumption  that  C  is  known  for  the  tender,  a  com¬ 
parison  is  made  between  this  value  and  C*.  If  C  >  C1  an  iterative 
solution  is  required  to  determine  the  X  which  results  in  C1  =  C. 

The  corresponding  z^’s  obtained  are  the  optimum  load  list. 

If,  on  the  other  hand,  C  ^  C  1  the  model  indicates  that  the  ten¬ 
der  cannot  effectively  meet  demands  without  re-supply  during  the 
period  under  consideration.  In  the  event  this  occurs  several  decision 
criteria  can  be  adopted.  If  demand  for  each  item  is  relatively  uni¬ 
form,  one  solution  would  consist  of  decreasing  the  tender  load  list 
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time  requirement.  For  example:  Let  C1  =  15,000  cu. ft. represent 
a  ninety-day  load  list  minimum  requirement.  If  C  =  5,000  cu.  ft.  , 
then  the  tender  can  in  effect  maintain  a  30  day  load  list  by  decreasing 
each  proportionally. 

This  criterion  has  the  advantage  of  maintaining  a  load  list  which 
consists  of  all  items  that  are  demanded.  The  most  obvious  disadvan¬ 
tage  is  the  increased  re- supply  requirement. 

Another  possibility  is  the  elimination  of  items  from  the  load 
list.  This  requires  a  decision  criterion  to  select  the  item(s)  which 

should  be  deleted.  Since  C1  decreases  as  \  increases  .the  model 

M  • 

indicates  that  the  item  corresponding  to  min(  )  should  be  the  first 

item  removed  from  stock.  If  (C1— c^  a^  =  C  ;  a^  <  z^)  then  an 

amount  (z^  -  a^)  can  be  maintained  for  this  item.  If(C!—  cjzi>C) 

then  the  item  corresponding  to  min  (  — ~  )  is  eliminated  completely. 

M  * 

A  new  min(  -gr-  )  is  selected  from  the  remaining  list  of  N-l  items 

and  the  process  is  repeated.  When  a  X  can  be  found  such  that  0=0* 

the  z^'s  have  been  determined  for  this  decision  criterion.  This 

particular  policy  insures  that  the  items  with  the  highest  essentiality 

to  cube  ratio  (__±_)  are  retained  on  board. 

C-L 

hhat  period  during  which  the  tender  must  be  able  to  operate 
without  benefit  of  re- supply  and  still  maintain  effective  support 
capability. 

2  Note  that  other  policies  could  have  been  chosen.  For  example, 
one  might  choose  a  policy  that  insures  that  items  with  high  values  of 
P  are  retained  on  board,  while  items  with  low  values  of  P  are  eliminated. 
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It  permits  tender  support  to  be  maintained  at  some  reduced  effective¬ 
ness  level.  (This  reduced  level  can  be  a  significant  one).  In  effect, 
this  is  analogous  to  an  unexpected  demand  situation  and  will  be  con¬ 
sidered  further  in  connection  with  tender  re-supply  requirements. 

The  iterative  nature  of  the  solution  and  the  number  of  items 
comprising  load  lists  indicates  the  necessity  for  utilizing  electronic 
computers  if  the  system  is  to  be  incorporated  effectively.  Appendix 
III  contains  a  computer  description  and  flow  chart  for  the  model. 
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5. 


Example  and  General  Results: 


A  numerical  example  will  serve  to  illustrate  the  procedure 
for  determining  an  optimal  tender  load  list  utilizing  the  proposed 
model.  Consider  the  case  where  the  demand  for  each  item  is 


characterized  by  the  normal  distribution.  From  equation  (5) 


subject  to  the  available  stowage  constraint.  It  should  be  noted  that 

each  item  could  have  different  demand  distributions  and  the  model 

would  still  be  valid.  The  only  requirement  for  utilizing  the  model 

is  that  the  density  functions  be  continuous.  * 

Table  (1)  contains  the  parameter  values  pertinent  to  the 

example.  In  accordance  with  the  procedures  of  the  previous  section, 

an  available  tender  capacity  of  10,152  is  assumed.  It  should  be 

apparent  that  due  to  the  iterative  nature  of  the  solution, 

will  be  approximately  equal  to  10,152  for  an  optimum  set  Z. 

(The  most  obvious  reason  is  the  need  for  rounding  off  the  values 

obtained  to  intergers.  In  order  to  insure  that  the  space  constraint 

is  met,  the  z^’s  obtained  should  be  rounded  off  to  the  lower  interger 

Ni  * 

with  the  result  that  ~  ). 

The  values  assigned  to  the  items  listed  in  Table  (1)  reflect 
varying  combinations  of  essentiality,  cube,  expected  demand  and 


^Another  model  could  be  developed  which  would  utilize  discrete 
functions . 
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standard  deviation. 


The  essentiality  valves  in  Table  1  were  arbitrarily  chosen  to 
show  that  items  1-8  are  approximately  four  times  as  essential  as 
items  9-16.  In  a  real  life  problem  involving  a  large  number  nf 
items,  a  scheme  using  the  military  essentiality  codes  could  be  adapted 
for  use  in  this  model. 


Table  1 


X 

T1 

h 

Etxa 

Hi 

NV, 

cc 

Ml 

P= 

<3!  -CL 

i  , 

100 

12 

? 

12 

.  75 

1.  33.  . 

6.25 

2 

100 

3 

? 

12 

.75 

1.33.  . 

25 

3 

50 

12 

9 

12 

.75 

1.33.  . 

3.  12 

4 

50 

3 

9 

12 

.75 

1.33.  . 

12.5 

5 

100 

12 

9 

3 

3 

.333. . 

Zb 

6 

100 

3 

9 

3 

3 

.333. . 

100 

7 

50 

12 

9 

3 

3 

. 333. . 

12.5 

8 

50 

3 

9 

3 

3 

.333. . 

50 

9 

100 

12 

2 

12' 

.  1666. . 

6 

1.  38 

10 

100 

3 

2 

12 

.  16^6.  . 

6 

- 

11 

50 

12 

2  ‘ 

12 

.  1666.  . 

6 

0.69 

12 

50 

3 

2 

12 

.  1666,  . 

6 

2.76 

13 

100 

12 

2 

3 

.666. . 

1.5 

5.55 

14 

100 

3 

2 

3 

.666. . 

1.5 

20 

15 

50 

12 

2 

3 

•  666. . 

1.5 

2.76 

16 

50 

3 

2 

3 

.666. . 

1.5 

11. 1 

In  order  to  determine  the  value  of  M  once  an  optimum  set  Z 
is  obtained,  equation  (3)  must  be  utilized.  For  the  case  when 
demand  is  characterized  by  the  normal  distribution, 

2-  oo 

(3)  StCza  =  P*i?l<^0&Xv  -V-  J  £u^cb<;  = 


12 


O') 


Tzw 


-^zvr 


Appendix  II  contains  a  detailed  derivation  of  equation  ( 3 1  - ) . 

Since  equation  (5)  requires  z.  0  this  implies  that  for  the 

case  of  the  normal  distribution  equation  (5)  and  ( 3 1 )  are  only  valid 
.  .  2 

when  - —  >  3.87.  (Fig.  1) 


Fig.  1:  Graph  of  the  normal  distribution  function 

Since  =>  Ziir^  >-3.8T. 

*  For  ‘Z.^>  O  ,  the  ratio  ^->  3  87  must  hold.  Examination  of 
Table  1  indicates  that  this  condition  has  been  satisfied. 

Due  to  the  requirement  above,  the  third  term  of  equation  (31) 


^  The  value  for-^-  is  dependent  on  the  accuracy  of  the  tables 
used.  The  value  of  3.87  was  obtained  from  the  Chemical  Rubber 
Company  Standard  Mathematical  Tables. 
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can  be  ignored  in  this  example  since  its  contribution  is  insignificant 
for  the  values  assumed. 


Table  2  is  a  presentation  of  the  results  obtained  when  the 
parameters  shown  in  Table  1  are  used  in  the  model.  The  optimal 
set  Zb=||^>,T-*.  ^  *  *  *^*\^  II  which  most  nearly  satisfies  the  arbitrary 
constraint  of  10,152  cu.  ft.  tender  capacity  was  found  by  using  A~.04\6o 
Note  that  items  9-12  generally  lag  in  percent  of  expected  average 
demand  filled  &ytaa,  for  each  value  of  A  .  This  is  intuitively 
appealing  because  these  are  items  with  low  military  essentiality  and 
high  cube.  The  items  with  high  essentiality  and  low  cube  (5-8)  are 
stocked  at  consistently  high  levels. 

Figures  2a  thru  2d  are  plots  of  percent  of  expected  average 
number  of  demands  filled,  versus  a  characteristic  number  P,  defined 


as 


Mi 

C  <Tc  ' 


Since  each  item  is  characterized  by  the  parameters  m,  <3“,  M  and  c, 
each  item  is  also  characterized  by  P.  Note  that  for  different  values 
of  C1  the  curve  shifts  somewhat  but  retains  the  same  general  shape. 
In  each  case  as  the  characteristic  number  increases,  the  percent  of 


expected  average  number  of  demands  filled  increases. 

The  effect  of  variability  of  demand  for  this  example  is  investi¬ 
gated  in  Figure  3  by  observing  the  behavior  of  items  9-12.  Since 
items  9  and  10  (also  11  and  12)  differ  only  in  their  values  of 
we  would  expect  any  differences  in  their  optimal  stock  levels  (  Z-C  ) 
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Table  2 


S-CC'Z.C  = 

8086 

A=  •  'ZSO 

Z6  -9341-  ^-MLS(ZiY«5l6 

T. 

2LL 

3(3*^ 

cc  Zo 

V\c  sCzCl 

■sU.^ 

<=-^1 

M.;s  h.i) 

1 

109.2 

98.46 

98.46 

1310.4 

886. 14 

111.6 

"C94 

98.94 

1339.2 

890.46 

2 

102.3 

99.61 

99.61 

1227.6 

896.49 

102.9 

99.73 

99.73 

1234.8 

897.57 

3 

59.  2 

48.46 

96.92 

710.4 

436. 14 

61.  6 

48;  94 

97.88 

739.  2 

440.46 

4 

52.3 

49.61 

99.22 

627.6 

446.49 

52.9 

49.73 

99.46 

634.8 

447.57 

5 

119.  1 

99.71 

99.71 

357.3 

897.39 

120.8 

99.79 

99.79 

362.4 

898.11 

6 

104.8 

99.92 

99.92 

314.4 

899.28 

105.2 

99.94 

99.94 

315.6 

899.46 

7 

69.  1 

49.71 

99.42 

207.3 

447.39 

70.8 

49.79 

99.58 

212.4 

448.11 

8 

54.8 

49.92 

99.84 

164.4 

449.28 

55.2 

49 . 94 

97.88 

165.6 

449.46 

9 

53.6 

53.55 

53.55 

643.2 

107.  10 

91.9 

90.11 

90.11 

1102.8 

180.22 

10 

88.4 

88.38 

88.38 

1060.8 

176. 76 

97.9 

97.52 

97.52 

1174. 8 

195.04 

11 

3  .  6 

3.55 

7.  10 

43.2 

7.  10 

41.9 

40.11 

80.22 

502.8 

80.22 

12 

38.4 

38.38 

76.76 

460.8 

76.76 

47.9 

47.52 

95.04 

574.8 

95.04 

13 

108.  1 

98.21 

98.21 

324.3 

196.42 

110.6 

98.76 

98.76 

331.8 

197.52 

14 

102.0 

99.55 

99.55 

306.0 

199.  10 

102.7 

99.69 

99.69 

308.  1 

199.38 

15 

58.  1 

48.21 

96.42 

174.  3 

96.42 

60 . 6 

48.76 

97.52 

181.8 

97.52 

16 

52.0 

49.55 

99.  10 

156.0 

99.  10 

52.7 

40.  ,9 

99.38 

158.  1 

99.38 

15 


2.cc.-z:  - 

> 

cnis- 

= .0833 

Z>\LsOzO) 

^bSS5 

^-cuz 

\  -  .O^-vb 

=  1015-2-  4582- 

0 

V. 

sUO 

^'U'sOz.O 

stzO 

1 

114.6 

99.35 

99.35 

1375.2 

894. 15 

119. 1 

99.71 

99.71 

1429. 2 

897. 39 

2 

103.6 

99.83 

99.83 

1243.2 

898.47 

104.7 

99.92 

99.92 

1256.4 

899.28 

3 

64.6 

49.35 

98.70 

775.2 

444. 15 

69. 1 

49.71 

99.42 

829.2 

447.39 

4 

53.6 

49.83 

99.66 

643.2 

448.47 

54.7 

49.92 

99.84 

656.4 

449.28 

5 

122.  9 

99.87 

99.87 

368.7 

898.83 

126.4 

99.94 

99.94 

379.2 

899.46 

6 

105.7 

99.96 

99.96 

317.  1 

899.64 

106.6 

99.98 

99.98 

319.8 

899.82 

7 

72.9 

49.87 

99.74 

218.7 

448. 83 j 

76.4 

;  49.94 

99.98 

229.2: 

449.46 

8 

55.7 

49.96 

99.92 

167.  1 

449. 64 | 

56.6 

49.98 

99.96 

169.8 i 

449.82 

9 

100.0 

95.21 

95.21 

1200.0 

190.421 

108.0 

98.21 

98.21 

1296.0  : 

196.42 

10 

100.0 

98.80 

98.80 

1200.0 

197.60 

102.0; 

99.55 

99.55 

1224.0 

199.  10 

11 

50.0 

45.21 

90.42  , 

600.0 

90.42 

58.0 

48.21 

96.42 

696.0 

96.42 

12 

50.0 

48.80 

97.60 

600.0 

97.60 

52.0 

49.55 

99. 10 

624.0 

99.  10 

13 

113.8 

99.25 

99.25  ‘ 

341.4 

198.50 

118.4 

99.67 

99.67 

355.2 

199.34 

14 

103.4 

99.81 

99.81 

310.2 

199.62 

104.6 

99.91 

99.91 

313.8 

199.82 

15 

63.8 

49.  25 

98.50 

191.4 

98.50 

68.4 

49.67 

99.34 

205.2 

99.34 

16 

53.4 

49.81 

99.62 

160.2 

99.62 

54.6 

49.91 

99.82 

163.8 

99.82 

16 


for  various  total  cube  requirements  to  be  due  to  variability  of  demand 

(cr^ . 

Figure  3  supports  the  argument  that  when  a  large  total  cube 
is  available,  items  with  high  values  of  <3"  will  be  stocked  in  larger 
amounts  than  items  with  small  <5~  .  However,  when  the  total  cube 
available  becomes  more  restricted,  the  stock  of  items  with  high 
variability  of  demand  is  reduced  rapidly.  Plots  for  other  items  which 
differ  only  in  their  value  of  G"  show  a  similiar  trend. 

Reference  to  Figure  4  indicates  that  as  total  available  cube  in¬ 
creases,  the  total  essentiality  of  the  load  list  increases  also,  but  at 
a  decreasing  rate.  Planners  who  wish  to  increase  the  spare  parts 
capacity  of  a  tender  would  have  to  weigh  the  cost  of  increasing 
capacity  against  the  value  of  the  added  cargo.  Use  of  a  computerized 
version  of  the  model  with  a  large  number  of  items  would  provide  a 
new  optimal  load  list  for  the  proposed  increased  capacity.  Reference 
to  a  curve  similiar  to  that  shown  in  Figure  4  would  be  an  aid  to  the 
planner  who  must  determine  whether  the  essentiality  of  the  added 
inventory  is  worth  the  cost  of  increasing  tender  spare  parts  capacity. 


17 


100 


o 

o 


o 

00 


o 

h- 


O 

vO 


O 

lT) 


Item  Characteristic  Number  (P) 
Fig  2(a) 


o 

o 


o 

00 


o 


o 

o 


Item  Characteristic  Number  (P) 
Fig  2(b) 


— DOT 


Item  Characteristic  Number  (P) 
Fig  2(c) 


o 

o 


o 

00 


o 

r- 


o 


o 

to 


o 


Item  Characteristic  Number  (P) 
Fig  2(d) 


50 


=1.00 

Kl.  -  Z. 

Xv.  =  12, 


75W 


shtt - sW 


L 


tcjWo- 


9000 9500 

TOTAL  CUBE 


10000 


Fig  3:  Effect  of  Variability  of  Demand 
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6.  POLARIS  Tender  Re -supply  (General) 

The  maximum  deterrent  value  of  the  FBM  weapon  system  can 
be  realized  only  if  the  maximum  number  of  SSB(N)'s  are  maintained 
in  the  designated  areas.  The  logistic  support  concept  is  based  on 
satisfying  this  requirement.  Lack  of  essential  spare  parts  can  be 
interpreted  as  a  deficiency  in  the  number  of  ready  missiles  scheduled 
to  be  on  station. 

A  load  list  is  developed  to  allow  the  tender  to  provide  effective 
support  for  a  given  period  of  time  in  the  event  re- supply  is  cut  off 
for  any  reason.  As  a  result,  replenishment  of  the  tender  is  in¬ 
fluenced  by  the  following  factors: 

(1)  maintain  stock  levels  at  or  near  the  optimum  load  list 
point ,  and 

(2)  provide  delivery  of  nonstandard  demands  (items  not 
included  on  the  load  list  for  one  reason  or  another) 
from  the  stock  point  to  the  tender. 

The  factors  considered  in  selecting  the  mode  of  transportation 
(air  or  surface)  for  an  item  include  the  transit  time  available  , 
priority,  physical  characteristics  of  the  material,  capability  and 
availability  of  the  various  carriers,  security  and  cost.  Of  these 
factors  those  necessary  to  effect  delivery  to  meet  program  require¬ 
ments  should  govern  and  considerations  of  cost  and  equitable  traffic 
distribution  should  be  regarded  as  secondary. 
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In  order  to  illustrate  the  considerations  mentioned  above  con¬ 


sider  the  case  where  a  tender  is  assigned  to  support  X  SSB(N)'s. 
The  load  list  for  the  tender  is  designed  to  provide  Y  days  support 
during  which  time  each  SSB(N)  is  alongside  for  Z  days.  (Fig  5). 

Y  can  be  interpreted  as  one  cycle  during  which  an  SSB(N)  has  an 
upkeep  period  of  Z  days  and  a  patrol  of  W  days  (Y  =  Z+W).  If  the 
assumption  is  made  that  the  alongside  availability  period  for  each 
SSB(N)  is  the  same,  then  3  K  is  the  number  of  SSBfN^s  at 

the  tender  at  any  time-  Under  this  assumption  an  SSB(N)  departs 
and  a  new  one  arrives  every  days.  When  demand  by  each 


SSB(N)  is  assumed  to  be  uniform,  stock  levels  aboard  the  tender 


where  there  is  no  re- supply,  the  tender  is  completely  depleted  of 
spare  parts  at  the  end  of  Y  days. 


Y 


Z.  — H-* 


-w 


H 


Tender  Support  Cycle 
Fig.  5 
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7. 


Routine  Re- supply 


Consider  the  first  factor  influencing  re- supply  mentioned 
previously.  The  average  percentage  of  the  optimum  amount  ('Z.'v.  ) 
which  the  tender  maintains  will  be  defined  as  the  effective  inventory 
level.  This  level  is  a  function  of  demand,  re-order  policy,  and  the 
lead  time  associated  with  the  transportation  of  the  item.  Consider 
the  case  where  demand  for  all  items  is  uniform  over  a  ninety-day 
period.  Fig  (6)  is  a  plot  of  effective  inventory  levels  at  the  tender 
versus  lead  time  required  for  delivery  as  a  function  of  re-order  policy. 

The  graph  can  be  utilized  in  several  different  ways.  If  a  re¬ 
order  policy  is  chosen  the  graph  indicates  the  inventory  level  that 
can  be  maintained  for  various  items  as  a  function  of  lead  time.  If 
lead  times  are  fixed  and  a  certain  inventory  level  is  desired  to  be 
maintained  then  the  re-order  policy  can  be  obtained  from  the  graph. 

It  should  be  noted  that  due  to  the  capacity  of  the  ship  and  the 
lead  time  required  for  delivery  of  an  item  the  tender  will  never 
average  an  amount  equal  to  the  load  list  quantity  for  that  item. 

Several  examples  will  serve  to  illustrate  the  use  of  Fig  (6). 

(a)  Consider  the  case  where  an  effective  inventory  level 
of  90  per  cent  is  desired  for  an  item.  The  graph 
indicates  that  this  level  can  be  maintained  by  re¬ 
ordering  when  the  stock  level  drops  to  90%  provided 
the  lead  time  is  10  days.  Another  possibility  is  to 
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re-order  when  stock  levels  drop  to  95%  if  lead  time 
is  16  days.  The  re-order  quantity  in  each  case  is  the 
demand  to  date  plus  expected  demand  over  the  lead 
time  . 

(b)  If  the  lead  time  for  delivery  of  an  item  is  fixed  then 

the  re-order  policy  determines  the  effectiveness  levels 
maintained.  With  a  lead  time  of  12  days  a  60%  policy 
results  in  an  effective  level  of  76%,  an  80%  policy 
results  in  a  level  of  84.5%,  a  90%  policy  in  an  89% 
level  and  a  95%  re-order  policy  in  a  level  of  93%. 

If  the  lead  time  is  30  days,  which  is  analogous  to  that  of  surface 
re- supply  by  AK  ^then  the  maximum  possible  effectiveness  level  is 
83%  for  the  example  shown. 

It  is  obvious  that  high  inventory  levels  can  be  maintained 
effectively  only  for  the  case  when  the  lead  times  are  short.  This 
illustrates  the  requirement  for  air  transportation  of  spare  parts  as 
opposed  to  surface  re-supply. 

In  addition  to  the  effect  on  effective  inventory  levels  maintained, 
air  re- supply  has  the  advantage  of  reducing  the  confusion  factor  at 
the  tender  when  supplies  are  received.  Rather  than  being  saturated 
with  parts  delivered  by  an  AK  all  at  one  time  the  tender  is  able  to 
maintain  better  control  over  spare  parts  storage.  This  in  turn  tends 
to  reduce  "loss"  of  items  due  to  misplacement. 
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Lead  Time 


8. 


Determination  of  resupply  method  for  NIS/Nonstandard  items. 
Consider  a  spare  part  that  is  required  by  an  SSB(N)  and  that  is 


not  in  stock  (NIS)  on  the  tender.  If  it  is  assumed  that  a  spare  part 
NIS  on  the  tender  is  always  available  at  the  stock  point  then  the 
mode  of  re- supply  transportation  (air  or  surface)  can  be  determined 
by  using  a  criterion  developed  as  follows: 


Let ,  - 

Julian  date  of  departure  of  AK  from  stock  point. 

= 

Days  AK  spends  in  transit  to  tender. 

li 

Days  remaining  until  j—  SSB(N)  starts  patrol. 

X) 

l#"- 

II 

Julian  day  that  j—  SSB(N)  generates  demand 
for  an  NIS  spare  part. 

c  = 

Days  required  for  loading,  unloading  and 
installation  of  spare  part. 

LTT  = 

Days  required  by  spare  part  to  travel  by  air 
from  stock  point  to  tender  (MATS,  commercial, 
etc). 

If  —  ^Av^-V-~^\vc  -V-  C*  ,  re- supply  the  spare  part  by  AK; 

If  Rj  -*-Sj  <  ^ak-V-Tak+C.  >  re-supply  the  spare  part  by  air, 
subject  to  the  constraint  '  • 

If  -v-S:  <  \_TT~  ,  then  it  is  impossible  to  re-supply  the  spare 

J  sj 

part  prior  to  the  SSB(N)'s  scheduled  departure. 
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9. 


Conclusions 


1.  The  number  of  stockouts  will  be  reduced  if  stocks  are 
maintained  in  accordance  with  the  model.  The  model  maximizes 
the  average  number  of  demands  filled  over  some  length  of  time. 
This  is  equivalent  to  minimizing  the  number  of  stockouts  over  the 
same  length  of  time. 

Z.  As  the  spare  parts  stowage  capacity  of  the  tender  is  in¬ 
creased,  the  total  worth  (military  essentiality)  of  the  inventory  is 
increased,  but  at  a  decreasing  rate.  For  example,  if  a  small 
tender  and  a  large  tender  each  have  their  cubic  capacity  increased 
by  the  same  constant  amount,  the  small  tender  has  a  larger  per¬ 
centage  gain  in  the  military  essentiality  of  cargo  than  the  large 
tender . 

3.  Items  with  high  essentiality  to  cube  ratios  are  generally 
carried  at  high  stock  levels.  For  those  items  with  nearly  the  same 
pattern  of  demand  the  model  maintains  high  inventories  of 

high  essentiality,  low  cube  items,  and  low  inventories  at  low 
essentiality,  high  cube  items. 

*  4.  When  items  must  be  deleted  from  inventory  because  of 

space  limitations,  those  with  the  lowest  essentiality  to  cube  ratio 
are  the  first  to  be  eliminated. 

5.  The  spare  parts  stowage  capacity  of  the  tender  is  the 
limiting  factor  in  determining  the  length  of  the  period  over  which 
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demands  filled  should  be  maximized.  The  load  list  should  be  designed 
to  maximize  demands  filled  over  a  given  length  of  time. 

6.  Re -supply  by  air  is  required  to  maintain  high  inventory 
levels  of  essential  spare  parts  aboard  the  tender.  In  order  to  main¬ 
tain  essential  spare  parts  inventories  close  to  the  optimum  load  list 
level  at  all  times,  it  is  necessary  to  have  short  re- supply  lead  times, 
or  a  large  tender  capacity.  Since  tender  capacity  is  limited,  high 
inventories  can  only  be  achieved  by  utilizing  the  short  lead  times 
associated  with  air  re- supply. 

7.  Re- supply  ships  cannot  meet  demands  for  non-standard 
items.  The  long  lead  time  associated  with  ship  re- supply  precludes 
the  possibility  of  meeting  demands  for  non-standard  items  as  they 
occur.  For  example,  it  is  a  rare  occurrence  when  a  submarine  can 
execute  an  order  for  a  non-standard  item  and  have  it  delivered  by 
ship  re-supply  before  going  on  patrol  again. 

8.  The  iterative  nature  of  the  model  and  the  large  number 

of  items  in  a  real  life  inventory  make  the  use  of  a  computer  mandatory. 
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appendix  I 


Derivation  of  Optimal  Tender  Load  List  Model 


(1) 

EXxJl  = 

J  Q  *w  d.x*v.  . 

(2) 

J,  U^--zl>)ftCXv.')cLx^  . 

^  V 

1 

(3) 

-  ’ 

)  *»-  jiCXvV&x;,  'z-i  j 

(4) 

ELLxa  - 

-v-  be  (z.0  . 

Essentiality  of  a  spare  part  is  an  important  factor  in  deciding 
on  a  load  list  which  will  maximize  the  average  number  of  demands 
filled  by  the  tender.  On  the  basis  of  past  usage  data  the  model 
should  indicate  the  quantity  of  each  item  to  be  stocked  such  that  the 
total  essentiality  function  M  is  a  maximum  subject  to  the  tender 
stowage  limitations. 

The  function  )  will  be  introduced  to  represent  a 

utility  valve  associated  with  the  ability  of  the  tender  to  meet  demand 
for  an  item. 


X  i  ^  ~Z.  I 
X  {_  >  ~2_  ^ 

L  -  \  j  •  •  •  *  Nl 


It  should  be  apparent  that  if  demand  is  less  than  what  is 
stocked  the  tender  is  able  to  provide  100%  service.  On  the  other 
hand  if  demand  exceeds  the  stock  level  then  the  tender  does  as  best 
as  it  can  by  providing  an  amount  ~Z-C.  • 
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In  accordance  with  the  measure  of  effectiveness  mentioned 


previously  the  objective  is  to  maximize  the  total  utility  value  as 
defined  above 

Max.  7. ......  Z.WV  ElLj- 

subject  to  the  constraint 

C.  ^"2  v  ^  C*  2-  ^  ■>  O  • 

^.Cx^-z^Y^  J - JCx*>yi..xN>|ci)(t...c\xKi 

where  \  y  X  i  . . is  the  joint  probability  function  of 

X^jX^^******Xki  • 

By  assumption  (7)  this  joint  probablity  density  function  can  be 
written  as: 

. x*\  «  . 

C<3 

LBv  definition  S9  fwUC>A-xL  =13 

As  a  result  can  be  expressed  as  follows: 

ECz^uc.xcQ  =  yo  . ■ 

Let  M;  S.fc't  = 

M;  ^  ^wCx;,2^fvCx^<Sx; 

Maximize  ^  v.  ^  v.C2. 

Subject  to:  -  C_ 

L 

The  problem  presented  is  one  of  maximizing  a  function  of  v\. 

variables  where  the  variables  are  related  by  a  single  constraint 
equation.  The  problem  could  possibly  be  solved  by  eliminating  some 
of  the  variables  using  the  constraint  equation.  Eventually  the  problem 
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could  probably  be  reduced  to  an  ordinary  maximization  problem. 
This  procedure  is  not  always  feasible  and  therefore  the  technique 
of  Lagrangian  multipliers  will  be  utilized. 

To  find  the  points  of  . where  the 

variables  "Z^Z*.. .  -"Z.^  are  related  by  the  equation  —  C=  O 

the  function  L  is  defined  as  follows: 


where  X  is  a  Lagrangian  multiplier 

=  —  Xcc  , 

c  v_=  \,  z_,. _ M 

As  a  result^the  following  set  of(n  +  Inequations  in(n  +  1^ 
variables  ‘Z.v>Zi...'ZNj  and  X  is  obtained. 


_ 


X  (c.»-z.v-  C_')  -  O 

o 


»4 


X  ■$—  (  cX  -  o 
-  C  -  o 


By  rearranging  the  terms  the  above  equations  can  be  put  in  the 


form 


(6) 


jeL_  \Lmc.  adz;.  VI  -  Xcc 
dz; 
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Recalling  equation  (3) 

(61) 

Lebnitz's  rule  is  utilized  to  evaluate  the  left  hand  side  of  this 
equation  since  the  units  of  integration  depend  on  the  parameter  'Zi  . 

Rule  :  V  WCXjVjdfc  =  >  tl  b'(V)  -  V\_Ca-CO , -\Qcx!co 

-4-  r  J^brlx.-'c*  cL* 

Jg-CM 

^_C5o  X.fvU^ckxw4-2c^  ^ 

=  C^w^u^-cm-o]  +  Co- 

=  L.  £l<^0<kx^ 

Since  r  frltxC) &-xi  ~  IL  -  X0  .fcCxO &xl  -t-  j*2,  ^^CxOcCxc. 

C*^5  "2.  * 

Then  J  £(jO&XC.  -  1  —  J“  iLOO&XC 

Z«L  o 

By  substituting  and  rearranging  terms  equation  (6 1  )  reduces 
to  the  following: 

1 —  SQ  5^bO<$b<v  - 

or  Jo  =  1-XC^.)  (5) 

Equation  (5),  subject  to  the  constraint  equation,  provides  the 
method  where-by  the  set  Z  can  be  determined.  The  two  equations 
determine  a  unique  set  of  *Z^'s  provided  that  for  all  L  ^  ^CX^>Q 
for  all  values  of  . 

The  distribution  function  ,  ,  of  a  numerical  valued 

random  phenomenon  is  defined  as  having  as  its  value,  at  any  real 
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number  “Z_ 


,  the  probability  that  an  observed  value  of  the  random 


phenomenon  will  be  less  than  or  equal  to  the  number  2L  . 

-  PLZ-  £  2L~^  .  Since  the  left  side  of  (5)  is  a  distribution  function 
the  following  condition  applies: 

o-  ■C'fct-'Oi'ioc  4  1  . 


In  order  to  solve  for  21 L  in  equation  (5)  we  must  assume  a 
value  for  A  .  It  would  materially  decrease  the  iterative  require¬ 
ments  if  the  possible  values  of  \  could  be  confined  to  a  portion 


of  the  real  line. 

Since 

O  S*'  £  1x0  Ax;.  & 

1 

Then 

°  <  i.-  X(  ^ 

>  X 

-X(^<  o 

>  X  ^  O 

• 

%  • 

>  A  >  o 

These  limits  apply  for  all  v_  and  therefore  A  must  fall  in 
the  range  .  O  in  order  for  all  items  to  be 

included  in  the  load  list. 
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APPENDIX  II 


Derivation  of  S^C ’Zl')  for  the  case  where  demand  follows  the  normal 


distribution: 


where 

oo 


but 


=  SQ  Xi£CX^Ax^^--2.i  ^  £oo<U;_ 

S0  £cxOAo<v  -  ~  S0  £(xMki  fcXxO 

.\  £  £u^xc  =  i  -  y*"  f^cxO^x; 

4v.  o 

y2'  S.U^ckxc  -  1-  XC^ 

■'*  5°~  J.GO&Xl  «0-~  C±-A(-^;)3} 

.%  ^Cz;^  =  y^xc£:cxo^; -V  -z-i  C\- a-  A-^oQ 

In  order  to  solve  for  an  explicit  expression  must  be 

found  for  i  *iSc  CjO&x;.  .  The  procedure  is  as  follows: 

■zt  z.  —  -1  '  * ^ ^ 

(ii- 1)  Sa  XvL^x^  <^-x:  ”  X,  ^  "  d-x: 


-i_  (  *y  —  y^jv 

\  —  Z-  X 

By  taking  the  derivative  of  -==.  S_  '- 


\7^V 


the 


following  expression  is  obtained: 


a  J_  ~ 

-^TV^" 


- -L-  -  (Xy-^wljk. 

^  ^  V  crv  rc^ 


=  —  l*v.-»Q 


_ L  /  Xv-  vv\v  \  x 

*  ^  “c^T  ) 


,  3-TT  ^ 

-  _  x:  ^  -4-  e_ 


Integrating  both  sides,  of  the  equation  between  the  limits  of  O 

4 

and  ~Z>  ^  : 

Wtw 

y  -* 


>.  /XC-wvCX^  z; 
-  x  v  <3-^  1  r 


Sx;^r 


x/x:-  K.V 
-  A  o-v.  >. 


-±VPJ 


r:  -1(^1 

J0  ~^TT^y  1 
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Rearranging  terms: 

x(xv-m-y 


-\xrr<5V 


d*. 


•z-c 

=  5, 


«•  ^ 


v>U 


-U 


c3o<  -Carr* 


Note  that  the  expression  on  the  left  is  the  same  as  equation 
II- 1  except  for  a  factor  OT^  .  In  order  to  obtain  the  same  ex¬ 
pression  multiply  both  sides  by  OT^  which  results  in 


s 


•A^ivCT; 


-TV 

o 


-te-Sr'V- 


X  #*e-  ^  -  X  A: 


<s\ 


r?-'-  ,  _  _xV—  » 

but  _  \  d  Sfe  e-  =  j=- 

J0  dx*  -Xirv  -varr 

-i-  f-zc-^OX* 

-  5=r  ' 


<S~c 

-\"wr 


tcx^lr 


—  z.  v  cr:-  * 


G- 


G~C  —  v  <rc 
'  Vi-TT 


Observe  that 


A.  (  Xv-^V~ 

_  V  <3"  ' 


but 


v  (X^-W\A  -ZC 

—  3L  V  CTv.  '  c  \  ^ 


Vv\  ‘  C~  '-Jb-  e,”  ^  ^  ^Xv  =  KvC\-X(~;^ 


Combining  the  above  derived  relationships  an  explicit  equation  for 
results : 


(3.)  vv\c.C'~  ^  x4ri''i 

.a;  -ttfcY  a-.-  ■ 

for  "Z.^  >0 
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APPENDIX  III 


A  computer  program  based  on  the  flow  diagram  illustrated  in 
Figure  7  was  written  in  the  Neliac  language  and  compiled  on  the 
United  States  Naval  Postgraduate  School  CDC  1604  digital  computer. 
The  Neliac  statements  are  found  in  this  appendix  following  Figure  7.  ^ 
Outputs  for  several  choices  of  parameter  values  confirmed  the  re¬ 
sults  shown  inS  ection  5. 

As  indicated  by  the  computer  flow  diagram  in  Figure  7  the 
first  step  in  the  operation  of  the  program  is  to  carry  out  the  routine 
input  procedures. 

The  next  step  is  the  computation  of  the  value  of  lambda  equal 

/  \ 

to  the  smallest  fraction,  v  Q_~r  )  .  Denote  this  minimum  ratio  as 

min  (,  .  This  value  of  lambda  is  used  to  compute  all  values 

of  AREA  =  (  ,i_ —  X  -pjr  )  for  i  =  1,  .  .  .N.  Each  ~Z.[_  is  computed 

from:  ~  x  fevf 


'VJ-vX.  ~ 


which  provides  an  approximation  of  the  normal  distribution  function. 
In  the  expression  above:  C  ~  ^>8 

where  "  \  -V  O-  ^  O . 

-V-  O.OZnH3cHS  (^-O.C003C\H%(^] 

-V  O.  * 


^Details  of  the  Neliac  language  may  be  obtained  from  Prof. 
R.  M.  Thatcher,  Operations  Research  Department,  U.S.  Naval 
Postgraduate  School,  Monterey,  California. 
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The  product  C^.'Z.C.  is  computed,  stored  and  summed.  The  result, 
<^L~Z-C  is  stored  in  the  cell  named  SUMCZ. 

This  sum,  ,  is  then  compared  with  the  contents  of 

Big  C  ,  the  specified  total  cube  available  for  the  storage  of  spare 
parts.  If  SUMCZ  >  Big  C,  then  the  cube  of  the  load  list,  for  the 
given  lambda,  is  greater  than  the  capacity  allowed  for  storage.  When 
this  occurs  the  item  with  the  lowest  essentiality  to  cube  ratio  is  dis¬ 
carded,  a  new  X  =  minl^c^  computed  from  the  remaining  items 
and  is  used  to  determine  a  new  SUMCZ.  This  new  value  of  SUMCZ 
is  compared  with  Big  C.  If  SUMCZ  remains  larger  than  Big  C  the 
item  with  the  lowest  essentiality  to  cube  ratio  is  discarded.  This 
procedure  is  repeated  until  SUMCZ  becomes  equal  to  or  less  than 
the  cube  available,  Big  C.  When  a  new  load  list  is  found  that  meets 
the  cube  constraint  then  the  essentiality  of  the  load  list  is  computed 
and  a  print -out  is  obtained  for  Z;  ^  A  ^  and 

A  table  of  the  items  discarded  is  also  provided. 

If  the  cube  of  the  load  list  computed  for  the  initial  value  of 
lambda  (  vvou^  )  is  less  than  the  capacity  allowed  for 

storage,  then  a  print -out  is  obtained  for  and 

In  order  to  find  the  load  list  that  most  nearly  meets  the  cube  constraint 
it  is  necessary  to  evaluate  for  several  assumed  values 

of  X 
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The  number  of  values  of  A  to  be  assumed  is  specified  in  the 
program  by  the  parameter  K.  If,  for  example,  K~20  then  the 
variables  and  will  be  evaluated  for 

W  V 

\  w  ' 

twenty  equally  spaced  values  of  A  between  min  ^j-r-  and  zero. 

In  this  case  twenty  print-outs  will  be  provided.  The  optimum  load  list 
is  taken  as  the  load  list  that  most  nearly  satisfies  the  cube  constraint. 
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Fig  1, 


Computer  Flow  Diagram 
42 


*0 
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Neliac  Statements  for  Computer  Program 


A 


A 


A 


A 


vo'sh-h  *  o 


B\g>  C.  — 
MX  C 


SvGMA(\^  - 


(00.0  *  O 
S’O.O  x  o 
Voo -O  *  o 
50.0  *  O 
\OQ.O  *o 
50-0  XO 
\oo.o  XO 
SO.O  XO 

v?_.o  x  o 

R.  O  x  o 
vn..o  x  o 

V2..0  X  O 
VL>0  X  O 
V2_.  O  *0 
V7L.  O  *  O 
V2..0  <0 


(QO.O  X  O  , 
S’O.O^O, 

(oo.o  x  o, 
so.  0X0, 

100.0  *  o  j 
50-0*0, 
VOO.  O  X  O  J 

coso.o  XO  , 

■2>  C>  X  O- 
3.0  X  o, 

3-0  x  o, 

■3.0  X  o  , 
■3.0  X  O  , 
0.0  X  o, 
3.0  XO, 
0303.0  XO  , 


VO  (U=A 


A.O  *  O  . 
Ro  X  O  * 
R.o  *  o  , 
RO  XO  , 

z.o  *  O  , 
z.o  *0, 

z.o  XO, 
T..O  *0  , 


°vo  XO 
<=\.0  X  O. 

RO  X  O  , 

°v  o  X  o. 

"5.-0  X  o, 

2-0  *  O, 

7-  O  X  O  . 

000Z..O  x  o  j 


v*-.o  X  o  , 

V3L.O  X  o  , 

3.0  x  o  3 

3.0  X  O  , 
V2_.0  X  O  . 
VT..O  X  O  ' 

(3.o  x  o  , 

3.0  XO  . 


VZ..O  XO  , 

VT..O  *  O, 

3.0  X  o , 
3.o  x  o, 
va..  o  x  o  , 

VL- O  *0! 
3. 0  *  O  , 
Oco'i.Q  X  O 


XTE>\  NRUtA  - 
nxM  -oxo 


^>^,^(0,-1,  S,3, 
\0,v\,vz.,  V3,  VA  ,  VS, 


\<K  -  20.  o  x  O 
XX-  000, 

LAV\BQ^  =  CCOO.OCcC  *  O  , 
sonc.x  -  oooooo .  0000  x  o  , 

SOKVVv/S^i.  —  OOCOOO.OOCO  X  O  , 

AP.OA  OM  =  O  X  O  J 

"2_  (v<=^  r  oooooo.  OC00  x  o  , 

S'ELCvC')  -  Oooooo  -  Oooo  X  O, 


A 

9 


Indicates  a  blank  card 
Indicates  a  card  of  all  nines 
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A 


A 

/\ 


A 


A 


A 

A 


A 


VJOC-  ~  0*0  ,  ^REAn^  =  Ooo.OoaCi^O, 

YXOGaTT  ~  —  .  S'  *  o  1 

K~  0*0,  Y  ~  o  *  O  ,  T"  =  ooo.  oooo  *  O  . 

svGv  4)  zj>r  —o*o 

E.KP1  ■=  O  *  O  , 

EXPZ  -0*0, 

SQRT^px  =  0*0 

TvJc|)pT  =  b^ssi^slonz.xo,  "XX^OQ, 

FAFv_AQ  J  K\G\  >  Mvvsi  , 

Fv_AS,  ^  TF\J\<Ei  j  XTt_k\  YnsCARO  AO|  ~  Ooo  , 

~rL  =0*0  .  tx.  -oxo, 

*4oux_  -  i~m  ~?  oooTrrn  ooo  &  , 

;  START”  i  ,  {;<>-,  \  , 

,  1  «  Pfctfi<5RAR\  t  _OArA\CtTb\  Pa^E-X  »  ,  ,  ,  i 
PRIMT-  PARWAj  ’ 

O — ^  VTAC*  — *-  TTVAG^  — »~  VAMv^r^  — K\ 

Sqft,TF  (  T\.>i<JiF»TL  •,  ^CSi'CtTZ.TTA  , 

F\KiO  Hiv4  '•A’-qbC.  ; 

o  —  *,  VA  , 

^£P€AT  ; 

TF  m~>A  MR  O'A^i  -  FAOV_\_  FA  -V  1.  — »•  v\  }  v\  =  Vvitvi  : 

_  Fnush  .  "R<E.FEfcr.  y 

\MOa1  /  c.  Cnj  - *-  LAMQCq  , 

O - 2-  VAXM  ''Uv5\  , 

F<jbR.  K  -  Otq  vAKi  —  i_ 

K.  IF  TCA  KiR  OG  -  MOU.  •.  Ov'SCA.Q.O  Vsitpc.  . 

vA’C»<T  /  £-C.*C“\  — a-  \m^c  <  oxm^oa  : 

V<  — *■  VAVMi  F\AG\  *  ) 

OvSCAFO  VA<£C.  > 


IF  MFLA.4  =  O  ■  UAVA-SSQA 


-*■  MIM 


MAA'Sk 


>  •) 


tOAPATvS.  ARSSA,  : 

O - -■>  5'JV\CHi  , 

F#>R  3T  -  Olx\  -  1_ 

(.  ‘XF  3ZF3A  MHCXl  -  Mi  Ou_ 


i-Av\QO A  *  C.CCX3  /xmCxCA  — ^  AftOA 

x — rA  -  cm  •  o—  S.<5>~? _ «.  -i—  1 


O-S 


~  Ars~A  'CxA 


1*.  O  • 

TF  A  CXT  =  O  *. 

T-  CXG  ^  O-S 

AROAt-xS  <-  o.ocsoojl  :  O  — *-A 
wormeX  (xrsam5  \TA 
O  —  T  — •*”  TX  C^MT. 

„  „  _  A  P.EACXl  -  O-  F  — ARtAti  3  , 

APvM-<^v\j5-  <_  o.  ooool  *.  O  — *»  'X,  cowrr.  ; 

HORHEX,  (.  AROA.VA5  ;  T  A  * 


OfcMTT. 


■AuOxh' 


O&VTT  I 

X  "V  \  '*n-*,**i5^  ^  ^ 

J  {.  «.  TX  \  -  V  ■>  XX  <  VJ  UMStlA  \  *  \ 

>  l_AvAfo\j/\  4  \j  A.QeA  \  -  V  >  &«eA  CXI  < 
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\J  AREALS  \ 

=  \  >  5  VJ  <~T  \  =  \  >1“  >  1 

KKCXl  -v-  SVGVAA  *  T~  - ^  *S-  CX3 

2.CX1  *  c.Cx.“X  -v  ^'o^e=a_  — *»  2vjvac_"Zl  , 
<£va\t  *.  ^  , 

,  ^■prvx'TT  «  Sut^ct  \  —  \  >  sonnc^sl  >  ^ 

32F  TFv_A<^  —  O  •-  1^.5T,  c.^MPuTS  - 


A 

A 


A 


A 

A 

A 


A 


VEST-;  ^xotAC_r2;  <-  B\GkC_  *•  C^nVuT^S^.  j 

1  — -  FVAQ  , 

K>G*  -V-  1.  - »>  5 

(A3TSA  F\_AGv  - =»r  \ — »  VTTXM  X>\SCAOJ^  X  ''■O.  , 

VO  -V  i_  — *.  VO  , 

NiUUV_  - *-  X\BA  VOVR  X.V_“i  j 

f\^0  VAXVsi  VOcfc>C_. 


a— ^  nr 

Q.^>VAPOT-<El 


's*fc  *.  ,  <>  ,  >  ?  >  , 

Ua\x  VJVJ  "2:  OvJ  S  »  ,  } 


O 


Sonwl  , 


) 


<Jjvo  :  f<£>r.  X  -  o(.i^  NiNi  —  i_  {_  xf  \TeitA  xoRtx.1  -  vauvao  5<\p, ; 

“SvC^mAvXxOl  /  SQRT  2.  FE  — S\Qc^2_PX  , 


VAK.C.X1  /  5  \<aRA  CX3  — ■—  X  X  x  - ®-  X  - 

RES^S  X  X  — ►  X  ,  BXPF  (X  EXP1), 

(I.CT3-  vaxcx! /s\gvaaCx3  — »■  Y* Y— »-Y, 
VA^PT-  S"  *  X  expr  (A  -,  expz'1 , 

m\cx3  *  AWCACXI  vz.cxn  *  U-o -Area tso') 

V  S\G,^2_PX  *  0<P1  —  S\5QZ.Pr  *  , 


A 


A 


A 


VO  Lr'i  *■  SUz  -V-  'S'OtAVOS'Z:  - ■*’■  J 

XI  A-  - Jta.  X  T~ 

L^xkTT  <  X.X  IW\  "2.CX1  U  >  i  ,  SN<\F-.>? 

X>\SCA«.0  TAT3UE  '- 

,  t 

.  ^  <-<•  xrev\  \  OxSdARO  \  TA^x-J£=_  ■>■?  >  5  J 

.  OC  u\\  X  VO  UN  VvJ 

1  v>u  e  ou\  m  vj  svgvaa 

^  i  , 

FtjbPO  X  -  OCX>VoG-l_ 

XVv^itA  "Ox  sciVUXO  ^XT\  •  X ^ 

IV\£XA  "\)V5<lAOXi  CXI  -V  X  — ^  orro-x  XASCAOJO,  CXI  , 
S  'PRxxO'T'  A  \  1TX.NA  XJx^c^avO  VV\  VO  \v\ 

CO  CtI  v\\  vvk  CT3  VU  5XGCAA 
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A 

pvnu'sva  ;  .  o 

A 

9 

5T 

TT  =  OO 


A 

9 


V>Rv\KTT  "PAVX.AWA  l  C. 

,  (  44  U\  X  VO  U  C  \J  ^  W\  S\<S^A  »,  > 

F<f)R  3“  -  oC  10  —  X. 

{_  3--V  i_  — *.  =T3“. 

L  P^kTT  A  V\  3T  U\  Ms]  CT3  \\  C  ^3^  \\\  H*  C3^  UV  CJl 

.. 


^>\  -  O  *  O 
*  0*0 


5 

A  \  ~  o  x  o 

■O  \  - , 

M  \  OM,  -  O  *•  o  J  \  Ok^(£va  -0X0, 
s  \  rtj.  -  o, 

\2^  \  <j><fa-czxT  -  o.  *  o  , 

R.9^  • 


^  AR'S  A 


IP  AvVUSA  =  O 


l. 

o- 


exrr^x 


XP  o  >  a«.v=a  o  A<i.eA  k  c>.S"  •.  owujfco. 
jn _ o  — ^ 

C|qt>^o^5.  :  “F /  s.Frrz-  — >-  A  ,  (A*,  ^ 

G%5“  X  \=LRJF"  —  AR4EA  - *•  Vviow\ 

A  *  A  A  ,  O  —  A  — A  , 

EXPF  (  A  •,  A)  . 

•r^cIdt^jpx:  x  A  — 

KiON\/T)tlM^M  — >  DO-TATK  , 

T"  —  XeOATX  — 

IDF  ~T  <C  O 


O 


y  t 


A 

A 

9 


D£\_taTW.  /I - r 

0-0.  00005  T5i  o.occ^s  *. 

\XxrrMX.  o^kpo-tts- , 
VS  •.  ,  l  XujssiAl.  \  \jaujs.  \ 

<£>F  \  AU-V=!A,  y  \Ki  \ 
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